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We study the electron spin relaxation in intrinsic and p-type (001) GaAs quantum wells by 
constructing and numerically solving the kinetic spin Bloch equations. All the relevant scatterings 
are explicitly included, especially the spin-flip electron-heavy hole exchange scattering which leads 
to the Bir-Aronov-Pikus spin relaxation. We show that, due to the neglection of the nonlinear terms 
in the electron-heavy hole exchange scattering in the Fermi-golden-rule approach, the spin relaxation 
due to the Bir-Aronov-Pikus mechanism is greatly exaggerated at moderately high electron density 
and low temperature in the literature. We compare the spin relaxation time due to the Bir-Aronov- 
Pikus mechanism with that due to the D'yakonov-Perel' mechanism which is also calculated from 
the kinetic spin Bloch equations with all the scatterings, especially the spin-conserving electron- 
electron and electron-heavy hole scatterings, included. We find that, in intrinsic quantum wells, the 
effect from the Bir-Aronov-Pikus mechanism is much smaller than that from the D'yakonov-Perel' 
mechanism at low temperature, and it is smaller by no more than one order of magnitude at high 
temperature. In p-type quantum wells, the spin relaxation due to the Bir-Aronov-Pikus mechanism 
is also much smaller than the one due to the D'yakonov-Perel' mechanism at low temperature and 
becomes comparable to each other at higher temperature when the hole density and the width of 
the quantum well are large enough. We claim that unlike in the bulk samples, the Bir-Aronov-Pikus 
mechanism hardly dominates the spin relaxation in two-dimensional samples. 

PACS numbers: 72.25.Rb, 71.10.-w, 67.57.Lm, 78.47.-fp 



I. INTRODUCTION 



Much attention has been given to semiconductor spin- 
tronics both theoretically and experimentally due to 
great prospect of the potential applicationsjii^ii^ The 
study of the spin relaxation/dephasing (R/D) in semi- 
conductors contains rich physics and is of great impor- 
tance for the device application. Three spin R/D mecha- 
nisms have long been proposed in Zinc-blend semiconduc- 
tors, i.e., the Elliott- Yafet (EY) mechanism)^ caused by 
the spin-flip electron-impurity scattering due to the spin- 
orbit coupling; the D'yakonov-Perel' (DP) mechanism^ 
which is due to the momentum-dependent spin splitting 
in crystal without a center of symmetry; and the Bir- 
Aronov-Pikus (BAP) mechanism'' which originates from 
the spin-flip electron-hole exchange interaction. Previ- 
ous researches have shown that, in bulk systems, the EY 
mechanism is important in narrow-band-gap and high 
impurity semiconductors; the DP mechanism is dominant 
in n-type semiconductors; and the BAP mechanism can 
have significant effect in p-doped semiconductors i^i^ It 
is known that, in heavily p-doped bulk samples, the BAP 
mechanism is dominant at low temperature whereas the 
DP mechanism is dominant at high temperature with the 
crossover temperature determined by the doping level. 
In bulk samples with low hole density, the BAP mecha- 
nism has been shown to be irrelevant lii^ii^ In addition, the 
hyperflne interaction induced spin relaxation is another 



possible mechanism. 

In contrast to the bulk systems, the relative impor- 
tance of the BAP and DP mechanisms for the electron 
spin R/D in two dimensional (2D) systems, especially in 
p-type 2D systems, is still not very clear, sometimes even 
confusing. In Ref. [ll|, extremely long spin relaxation 
time (SRT), which is two orders of magnitude longer 
than that in the bulk sample with corresponding accep- 
tor concentrations, was reported by Wagner et al. in 
p-type GaAs quantum wells (QWs). The authors argued 
that the BAP mechanism is dominant at low tempera- 
ture. However in Ref. the SRT in p-type QWs was 
reported to be a factor of 4 shorter than that in compara- 
bly bulk GaAs by Damen et al. at low temperature. The 
authors also referred the BAP mechanism as a cause for 
the decrease of the SRT. Hence, two opposite experimen- 
tal results arrive at the same conclusions regarding the 
importance of the BAP mechanism. Moreover, Gotoh et 
al. further pointed out that the BAP mechanism should 
not be ignored even at room temperature.^^ They investi- 
gated the electric field dependence of SRT and found that 
the SRT decreases with the increase of the bias. They 
concluded that the decrease is from the BAP mechanism 
as the SRT due to the DP mechanism does not change 
with electric field. Actually, they overlooked the fact that 
the Rashba spin-orbit coupling^^ can also lead to the spin 
R/D due to the DP mechanism. Therefore, we believe 
that the decrease of SRT in their experiment cannot be 
a proof of the importance of the BAP mechanism. Very 
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recently it was shown that the SRT at room temperature 
can be increased at the (100) GaAs surface due to the 
relatively lower concentration of holes at the surface and 
the mechanism for the SRT was referred to as the BAP 
mechanism. Theoretically, Maialle^*' pointed out that 
the effect of the BAP mechanism in 2D systems is a little 
smaller than that of the DP mechanism at zero tempera- 
ture by using the Fermi golden rule to calculate the SRT 
in which the elastic scattering approximation was applied 
and consequently the nonlinear terms of the electron-hole 
Coulomb scattering were neglected. The SRT due to the 
DP mechanism (top) was also calculated by using the sin- 
gle particle approachfi*^ The author compared r^p and 
Tbap for different electron momentums (kinetic energies), 
and showed that these two SRTs have nearly the same 
order of magnitude in heavily doped QWs. However, 
Top calculated in Refs. 16] is quite cursory, because, un- 
der the framework of single particle theory the carrier- 
carrier Coulomb scattering, which is very important to 
spin R/Dfilii^ii^iS£i^ is not included. Also the counter 
effect of the scattering to the spin R/D is also not fully 
accounted ii^ii^i2£ii22i2^ Moreover, it is also important to 
calculate the spin-flip electron-hole exchange scattering 
explicitly in order to find out the effect of the nonlinear 
terms ignored in the Fermi golden rule approach by Ma- 
ialle et 26^27 -y^g g^jg^ want to find out the tempera- 
ture dependence of the relative importance of both mech- 
anisms in 2D systems, which to the best of our knowledge 
is still absent in the literature.. 

In order to accurately investigate the relative im- 
portance of the DP and the BAP mechanisms beyond 
the single-particle Fermi golden rule approach, we use 
the fully microscopic approach established by Wu et 
al^ by constructing and numerically solving the kinetic 
spin Bloch equations. In this approach, 
all the corresponding scatterings such as the electron- 
acoustic (AC) phonon, electron-longitudinal optical (LO) 
phonon, electron-nonmagnetic impurity, and electron- 
electron Coulomb scatterings are explicitly included. The 
results/predictions obtained from this approach are in 
very good agreement with varies experiments ;^°i^^'^^i'^° 
It was previously pointed out that, in the presence of 
inhomogeneous broadening, any type of scattering, in- 
cluding the Coulomb scattering, can give rise to the spin 
R/D.iiii^ii^i2ai2^ In this paper, in addition to the all the 
above mentioned scatterings in rt-type QWs as consid- 
ered in Ref. , we further add the spin-conserving and 
spin- flip electron- heavy hole Coulomb scatterings, both 
contributing to the DP mechanism and the latter further 
leading to the spin R/D due to the BAP mechanism. By 
solving the kinetic spin Bloch equations self-consistently, 
we obtain the SRT from the BAP mechanism from a fully 
microscopic fashion. We further investigate the relative 
importance of the BAP and DP mechanisms in 2D sys- 
tems. 

This paper is organized as follows. In Sec. II, we 
construct the kinetic spin Bloch equations and present 
the scattering terms from the spin-conserving and spin- 



flip electron-hole Coulomb scatterings. We also discuss 
the SRTs due to the BAP mechanism from different ap- 
proaches. Then we present our numerical results in Sec. 
III. We study the SRT due to both the DP and the BAP 
mechanisms under various conditions such as temper- 
atures, electron/hole densities, impurity densities, and 
well widths. We conclude in Sec. VI. 



II. KINETIC SPIN BLOCH EQUATIONS 

We construct the kinetic spin Bloch equations in intrin- 
sic and p-type (001) GaAs QWs by using the nonequilib- 
rium Green's function method:-^ 



Pk,(T(T' — Pk.o-cr' Igoh Pk,(T(T' Iscatt 
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with pk^air' representing the single particle density ma- 
trix elements. The diagonal and off-diagonal elements 
of Pk.iTo-' give the electron distribution functions /ko- and 
the spin coherence pk^cr-o-, respectively. The coherent 
terms Pfe.crcr' Icoh describe the precession of the electron 
spin due to the effective magnetic field from the Dres- 
selhaus ternt^^ r2(k) and the Hartree-Fock Coulomb in- 
teraction. The expression of the coherent terms can be 
found in Appendix A (and also Ref. [iB|). The Dressel- 
haus term can be written asi^ 
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in which (fc^) represents the average of the opera- 
tor — {d/dzY' over the electronic state of the low- 
est subband,-*' and 7 is the spin splitting parameter— 
which is chosen to be 11.4 eV-A'^ all through the 
paper. Pfe.o-o-' Iscatt ™ ^Iq. (H]) denote the electron-LO- 
phonon, electron- AC-phonon, electron-nonmagnetic im- 
purity, and the electron-electron Coulomb scatterings 
whose expressions are given in detail in Appendix A 
(see also Refs. fislfiol ioj ) . All these scattering are calcu- 
lated explicitly without any relaxation time approxima- 
tion. Moreover, we further include the spin-conserving 
and spin-flip electron-heavy hole scatterings as what fol- 
lows. 

The Halmitonian of electron-heavy hole interaction is 
given by 



H. 



eh = E 

k,k',q,cr = ±l,cr' = ±l 



e'»,9Ck+q,5Ck,f 6j 



k'-q, 
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where c (c^) and b (6^) are the annihilation (creation) 
operators of electrons in conduction (heavy-hole valence) 
band respectively. We denote a {a') to be ±1 throughout 
the paper. The screend Coulomb potential under the 
random-phase approximation reads^ 



eh^q 



e(q) 



(6) 
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with the bare Coulomb potential vq = 47re^/(5^ and 
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is the electron-hole plasma screening.-^^ In these equa- 
tions = + ql and /^^ (/k.tr) denotes the heavy 
hole (electron) distribution function with spin |cr (|cr). 
The form factors can be written as: 

fe{q.) = J dzdz'U^)Uz')e''''^'-''\c{z')Uz) , (8) 

h{qz) = j dzdz'rjh{z)rjh{z')e''i'^'-''')rj,,{z')rj,,{z) , (9) 



where ^c{z) {r]h{z)) is the envelope function of the 
electron (heavy hole) along the growth direction z.— 
The scattering term of this spin-conserving electron-hole 
Coulomb scattering can be written as: 
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where Pk = Pk i-i = Pl- _i i ■ This spin-conserving scat- 
's 2 2 2 

tering only enhances the total scattering strength mod- 
erately and contributes to the spin R/D due to the DP 
mechanism. 

The Hamiltonian of the spin-flip electron-heavy hole 
exchange interaction reads 



where AE^t is the longitudinal-transverse splitting in 
bulk, |(/'3Z5(0)p = l/(7raQ) is the 3D exciton state at 
zero relative distance, and — + icrky. For GaAs, 
AElt = 0.08 meV and ao = 146.1 A respectively,^ The 
form factor can be written as: 
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The matrix elements in the Hamiltonian are given by^ 
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(14) The scattering term from this Hamiltonian reads 
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If we denote K = k + k' as the center-of-mass momen- 
tum of the electron-hole pair, the product of the matrix 
elements in Eqs. (IT51) and (fTT)) can be reduced to: 
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It is noted that the spin R/D of the photo-excited holes 
is very fast^'^ and the electron-hole recombination is very 
slow compared to the electron spin R/D. Therefore, we 
take the hole distribution in equilibrium Fermi distribu- 
tion and /^^ = = /k- Further, by subtracting 
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in Eq. p6|) . one obtains: 
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In above equation, the terms A/k[(l — /k')/k'-q + |(/k' ^ 
/k'-q)(/k-q,+i) + /k-q,-i)] dcscrlbc the forward scatter- 
ing and correspondingly the terms A/k-q[/k/ (1— /k'_q) — 

5(/k' -/k'-q)(/k,-(-i + /k,-i)] describe the backward scat- 
tering. The SRT due to the BAP mechanism from the 
Fermi golden rulc^^ can be recovered from Eq. by ap- 
plying the elastic scattering approximation: e 
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Under this approximation, the nonlin- 



ear terms (in the sense of the electron distribution func- 



tion) iA/k(/k^ -/k^_J(/k- 
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2"JKWk' J]i'-qJ(Jk~ci, + l -r Jk-q,- 

ward scattering and iA/k-q(/k/ ^ /k'-q)(/k,+i + /k -i) 
in the backward scattering tend to zero. In the remain- 
ing linear terms, A/k = -^{<l>i/2 - 4>-i/2) ~ A/k_q 
with /ok = ^/.(e,,-^)^! by choosing /k,^ = 
Therefore, one recovers the SRT due to the BAP mech- 
anism from the Fermi golden rule approach>i^ 



1 
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In the next section, we will discuss the applicability of 
above equation which relies on the elastic scattering ap- 
proximation. 

In this work, we do not use the SRTs from the single- 
particle approach for both the BAP and DP mechanisms. 



Instead, we solve the kinetic spin Blocli equations self- 
consistently with all the scattering explicitly included. 
The detail of the numerical scheme is given in Refs. [l^, 
[20j . The spin relaxation and dephasing times can be 
obtained from the temporal evolutions of the electron 
distribution functions /k,cr and the spin coherence /Ok.iT-o- 
respectivelyi^i^ We will show that the SRT due to the 
BAP mechanism obtained from the kinetic spin Bloch 
approach can give markedly different results compared 
to the one calculated from Eq. ([20)1 by using the elastic 
scattering approximation, similar to the situation of the 
SRT due to the DP mechanism which has been discussed 
in great detail in our previous workS] ^^'^^i^° 



III. NUMERICAL RESULTS AND ANALYSIS 

The SRTs calculated from the kinetic spin Bloch equa- 
tions are plotted in Figs. 1 to 6. In these figures, the solid 
curves represent the SRTs due to the BAP mechanism 
(tbap) which are calculated from the kinetic spin Bloch 
equations by setting the DP term r2(k) = 0; the dashed 
curves are the SRTs due to the DP mechanism (r^p) 
which are calculated by setting 9/9k,o-CT'/^^lBAP — 0; 
the dash-dotted curves represent the total SRTs (rt„t^i) 
obtained from Eq. ([1]) with all the terms explicitly in- 
cluded. We always use different color and width of curves 
for different conditions. 
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FIG. 1: (a) SRT due to the BAP (solid curve) and DP (dashed 
curve) mechanisms and the total SRT (dash-dotted curve) vs. 
temperature T in intrinsic QW when a = 20 nm, electron and 
hole densities n — p = 2no, and impurity density rii = n. 
(b) (color online) SRT due to the BAP mechanism with full 
spin-flip electron-hole exchange scattering (solid curves) and 
with only the linear terms in the spin-flip electron-hole ex- 
change scattering (dotted curves) at different electron densi- 
ties against temperature T. no = fO^'^ cm~^. 



We first discuss the SRT in an intrinsic GaAs QW 
confined by Alo.4Gao.6As barriers. In Fig. [IJa), we plot 
the temperature dependence of the SRT for a QW with 
well width a = 20 nm. The electron (heavy hole) den- 
sity n (p) is 2 X 10^^ cm~^ and the impurity density 
rii = n. It is seen from the figure that the SRT due to 
the BAP mechanism is much larger than that due to the 
DP mechanism. Moreover, Tbap decreases dramatically 
with T at low temperature, followed by a more moderate 
decrease at high temperature. The temperature depen- 
dence of Tbap Can be understood as follows. When the 
temperature increases, more electrons and holes tend to- 
wards the lager momentum, hence the larger center-of- 
mass momentum K. This leads to a larger the matrix 
element in Eq. IjlSp . and consequently a larger scattering 
rate. Furthermore, the Pauli blocking which suppresses 
the scattering decreases with the increase of tempera- 
ture. Both leads to the decrease of the SRT due to the 
BAP mechanism. The temperature dependence of the 



SRT due to the DP mechanism has been well discussed 
in Refs. [UllillOIlll . Therefore we will not discuss the 
DP mechanism in detail in this paper. 
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FIG. 2: (color online) SRT due to the BAP (soUd curves) and 
DP (dashed curves) mechanisms and the total SRT (dash- 
dotted curves) vs. temperature T in intrinsic QWs at different 
densities (n = p = 2, 4, 6no) when a = 20 nm and rii ~ n. 
no = 10" cm"^. 

In order to see the difference of the SRT due to the 
BAP mechanism calculated from the full spin-flip scat- 
tering [Eq. (|17p ] and the one from the Fermi golden rule 
[Eq. (1^ ]. i.e., neglecting the nonfinear terms in Eq. (fT7)l . 
we plot the BAP SRT calculated from the Bloch equa- 
tions with only the linear terms in the spin-flip scattering 
as dotted curves for two different electron (hole) densities 
in Fig.[TJb). It is noted that for high electron density, the 
SRT due to the BAP mechanism from the Fermi golden 
rule is much smaller than Tbap at low temperature. Fur- 
thermore, the lower the temperature and/or the larger 
the electron density, the larger the difference is due to 
the "breakdown!^ of the elastic scattering approximation 
at low temperature and/or high density. This is in good 
agreement with the condition for the elastic scattering. 
The difference can be very small when the electron den- 
sity is smaller than 5 x 10^° cm^^ according our calcu- 
lation. Consequently the SRT for high electron density 
obtained in Ref. at zero temperature is much smaller 
than the actual one. Therefore, the effect of the BAP 
mechanism for high electron density at very low temper- 
ature is smaller than that claimed by Maialle et al. In 
fact, it can even be ignored. We further stress that the ef- 
fect of the BAP mechanism at low temperature and high 
electron density is far exaggerated in the literature due 
to the neglection of the nonlinear terms in the spin-flip 
electron-hole exchange scattering. 

In addition, in the presence of inhomogeneous broaden- 
ing, any scattering can give rise to spin R/D i-'^^'^^i^^i^'^i^'^ 
It is intuitive that the SRTs should satisfy: 



'^total '^DP "^BAP '^DP "^BAP '^differ 
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where Tbap is directly caused by the spin-flip electron- 
hole exchange interaction, r^p is from the inhomogeneous 
broadening when there is no spin-flip electron-hole ex- 
change interaction, and Tj^p corresponds to case with the 
presence of the spin-flip electron-hole exchange scatter- 
ing. The difference between and —J— is noted as 

^ TDP Top 

— - — . In our calculation we found — - — is so small 

Tj^iffcr differ 

that can be totally ignored. This is because the spin-flip 
electron-hole scattering is much smaller than the other 
scatterings. 
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FIG. 3: (color online) SRT due to the BAP (solid curves) and 
DP (dashed curves) mechanisms and the total SRT (dash- 
dotted curves) vs. temperature T in intrinsic QWs for dif- 
ferent well widths (a — 10 and 20 nm). n = p = 2no and 
impurity densities {rii — 0.5n and n). Note that the solid 
curves with the same well width but different impurity densi- 
ties exactly coincides with each other, no = 10^^ cm~^. 

Then, we discuss the temperature dependence for dif- 
ferent electron densities in intrinsic QWs in Fig. [21 One 
can see that Tbap decreases with increasing densites at 
high temperature but it behaves oppositely at low tem- 
perature. On the other hand, t^p decreases with increas- 
ing densities at all temperatures. We again interpret the 
density dependence of BAP mechanism by using the pre- 
vious arguments: at low temperature regime, i.e., in the 
degenerate limit, the Pauli blocking is enhanced by in- 
creasing the carrier density and/or lowering the temper- 
ature. Therefore, the scattering can be suppressed by 
increasing density. This causes an increase of Tbap- At 
high temperature regime, i.e. in the nondegenerate case, 
higher momentum states are occupied for larger density. 
This leads to a stronger scattering and hence Tqap de- 
creases with electron density. From this, we find that 
the relative importance of the DP and the BAP mecha- 
nisms does not change so much by changing the electron 
density. 

In Fig. [21 we plot the temperature dependence of the 
SRTs in intrinsic QWs for different impurity densities 
and well widths. It is clear that Tbap does not depend 
on impurity density, in other words, the curves corre- 
sponding to different impurities concentrations exactly 



coincide. However, Tdp can be enhanced due to the in- 
creased impurity scattering strength. If we enlarge the 
well width, both Top and Tbap become larger. This is be- 
cause the electron- hole exchange strength is weakened by 
the form factor Eq. (|15[) in the scattering matrix elements 
in the BAP mechanism for wider QWs. The leading term 
(linear term) of the Dresselhause spin-orbit coupling in 
Eqs. ([2J3]) is smaller for wider QWs in the DP mecha- 
nism. The variation of Tqp is larger than Tbap that is to 
say the relative influence of the BAP mechanism becomes 
more important for wider QWs. 
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FIG. 4: SRT due to the BAP (solid curve) and DP (dashed 
curve) mechanisms and the total SRT (dash-dotted curve) vs. 
temperature T in p-type QW when a = 20 nm, n = 0.5no, 
po = 4no, and rii — n. no = 10^^ cm^^. 

From our detailed investigations, we conclude that 
Tbap in intrinsic GaAs QWs is always larger than Top. 
At very low temperatures, the BAP mechanism can be 
ignored. However, it should be considered at higher tem- 
peratures for accurate calculating. Moreover, the relative 
importance of the BAP mechanism is increased by raising 
the impurity density and the well width. 

We now turn to study the SRT in p-type QWs. In Fig. 
[H we choose the well width a = 20 nm, n — 0.5 x 10^^ 
cm~^, p = n + Po = n + 4: X 10^^ cm~^, and rii = n. 
One can see that the magnitudes of t^p and Tbap are 
very close around T — 150 K. In p-type QWs, both the 
spin-conserving and spin-flip electron-hole scatterings are 
greatly enhenced by increasing the hole density. The for- 
mer gives rise to the increase of Tdp in the strong scatter- 
ing limi t^^i^^ and the latter gives rise to the decrease of 
Tbap- Therefore both SRTs are getting closer for larger 
hole concentration. In the case of Fig. [H the contribu- 
tions from the DP and BAP mechanisms are nearly the 
same around 150 K, and at lower and higher temper- 
atures, the contribution from the DP mechanism is no 
more than one order of magnitude larger than the BAP 
one. In addition, l/raiffor is still very small and can be 
totally ignored. 

We now analyze the temperature dependence of the 
SRT for different electron and hole densities in p-type 
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FIG. 5: (color online) SRT due to the BAP (solid curves) and 
DP (dashed curves) mechanisms and the total SRT (dash- 
dotted curves) vs. temperature T in p-type QWs with a = 20 
nm at different electron densities (n — 0.5 and Ino) and hole 
densities po = 2 and 4no. rii = n. no = W^^ cm~^. 
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FIG. 6: (color online) SRT due to the BAP (solid curves) and 
DP (dashed curves) mechanisms and the total SRT (dash- 
dotted curves) vs. temperature T in p-type QWs with n — 
0.5no, Po = 4no at different well widths (a — 10 and 20 nm) 
and impurity densities (ui — n and 2n). no = 10^^ cm~^. 

QWs. In Fig.O the calculated SRT for different electron 
and hole densities are shown. In Fig.[Bl a similar analysis 
is made for different well widths and impurity densities. 
The general features can be understood from the follow- 
ing. When the electron density becomes larger, both Top 
and Tbap become smaller with similar amplitude. (Note 
that n — 0.5no and uq are both within the nondegen- 
erate limit.) When hole density gets larger, both r^p 
and Tbap become smaller with the amplitude of the lat- 
ter being larger than the former (i.e., the importance of 
the BAP mechanism gets increased). This is because 
the electron-heavy hole scattering is markedly enhanced 
with the hole density. As the BAP mechanism is de- 
termined by the hole density, Tqap is very sensitive to 
the hole density. Nevertheless, Top is less sensitive as it 



is also determined by all the other scatterings. When 
the well width gets larger, r^p is enhanced with a larger 
amplitude at low temperature and with a small ampli- 
tude at high temperature, whereas Tqap becomes larger 
moderately. These results are similar to Fig. [3l Conse- 
quently the BAP mechanism becomes more important, 
especially around T = 150 K in the present case. When 
the impurity density gets larger, r^p becomes larger and 
Tbap does not change. This makes the relative effect of 
BAP mechanism become larger. 

From above features, we emphasize that the BAP 
mechanism is important in p-type QWs, especially for 
large well width and/or large hole densities (i.e. heavily 
doped) and large impurity densities. It is very differ- 
ent from the bulk systems in which the BAP mechanism 
is absolutely dominant at low temperature. Therefore, 
both the BAP and the DP mechanisms should be con- 
sidered to get the right SRT in QWs. 



IV. SUMMARY 

In summary, we have investigated the SRT due to 
both the DP and BAP mechanisms in intrinsic and p- 
type GaAs (001) QWs by constructing and numerically 
solving the fully microscopic kinetic spin Bloch equa- 
tions. We consider all the relevant scatterings such as 
the electron-AC phonon, electron-LO phonon, electron- 
nonmagnetic impurity, and electron-electron Coulomb 
scattering. Furthermore, the spin-conserving electron- 
heavy hole scattering, which enhances the total scatter- 
ing strength and therefore Tdp , and the spin- flip electron- 
hole exchange scattering, which induces the BAP SRT, 
are also included. 

We stress it is very important to calculate the SRT 
from our fully microscopic approach, especially at high 
electron density and low temperatures where the non- 
linear terms in the electron-hole exchange scattering be- 
comes very important. The SRT obtained from our fully 
microscopic approach is much larger than that from the 
Fermi golden rule. This means that the BAP mechanism 
is negligible at very low temperature and high electron 
density. We speculate this is also true in the bulk case. 
This is very different from the predictions in the litera- 
ture. 

We investigate the temperature dependence of the 
SRTs: The SRT due to the BAP mechanism Tbap de- 
creases rapidly with increasing temperature at very low 
temperature and slowly at higher temperature for both 
intrinsic and p-type QWs. It also decreases with electron 
density for both intrinsic and p-type QWs. For p-type 
semiconductors, it further decreases with hole density. 
We also compare the relative importance of the SRTs 
from the BAP and DP mechanisms. The SRT from the 
DP mechanism is also calculated from the kinetic spin 
Bloch equations which give the SRT also quite different 
from that from the single-particle approach as discussed 
extensively in our previous works J^'iSi'^iSii^ We find in 
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intrinsic QWs. the effect of the BAP mechanism is much 
smaller than that from the DP mechanism at low temper- 
ature and it is smaller by nearly one order of magnitude 
at higher temperature; In p-type QWs, the SRT from 
the BAP mechanism is comparable with the one from 
the DP mechanism around certain temperature (such as 
150 K in the case we study), especially when the hole 
density and/or the width of the QWs are large. For both 
the intrinsic and p-type QWs, the contribution from the 
BAP mechanism at very low temperature are negligible. 
We conclude that the spin R/D in QWs is very different 
from the bulk samples. In 2D case the BAP mechanism 
hardly dominates the spin relaxation. Instead, it is either 
smaller or comparable to the DP mechanism. 
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APPENDIX A: COHERENT AND 
SPIN-CONSERVING SCATTERING TERMS IN 
KINETIC SPIN BLOCH EQUATIONS 



The coherent terms can be written as 



dt 
dt 



= —a 



cob 



coh 



^^a;(k)Impk + fiy (k)Re/9k + 2alm ^ Vee,gPl+qPk 



(Al) 



jria;(k) -f Oj^(k) {fk, + l ~ fk-l) + i^Vee,q {h+q, + l ~ fk+q-l) PU - Pk+qifk, + l - fk-l) (A2) 



where V^^„ = 



The electron-impurity scattering terms read 

9fk,CT 

im q 

= {'^"^ ^9 '^(^k - 4-q) [(/k,-M + /k,-l)/Ok-q - (2 - fk-q,+l " /k-q,-l)pk] } - {k ^ k - q} (A4) 



dt 



[-2Trni C/2<5(4 - [/k,.(l - /k-q,.) - Re(pkP^_q)] } - {k ^ k - q} 



(A3) 



dpk 
dt 



in which |k ^ k - q| stands for the same terms previ- be 1 in our calculation) the charge number of the impu- 

1- ^ ' rity. The electron-phonon scattering terms are 

ously in s > but interchanging k ^ k — q. In these equa- 
tions C/q = I]„,(^i«Q/e(q))^/e(a0) with Zi (assumed to 



9fk,a 



dt 



dpk 

dt 



ph 



ph 



{-277 ^ gl,^,xS{el - 4_q - J^q,.,A)[^q,.,A(/k,a - /k-q,.) + /k,.(l - /k-q,a) 

-Re(pkPk-q)]} - {k ^ k - q| , 

H - 4-q - ^q«„A)[Pk-q(/k,-M + /k,-l) + (/k-q,-M + /k-q,-l - 2)pk 

qiJj.A 

-2A^qg,,A(Pk - Pk-q)]} - |k ^ k - q| , 



(A5) 



(A6) 
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where A represents the phonon mode. For the 
electron-longitudinal-optic-phonon (LO) scattering, the 
matrix element ff^ lq = /[[q^ + - 

Kq ^)/e((72); for clcctron-acoustic-phonon scattering due 
to the deformation potential, 5Q,de/ = ll^^^/el'S'z); 
and for that due to the piezoelectric coupling, = 
327r he e^^ '^^fJ^J^^I^^ feiQz) for the longitudinal phonon 



and 



9Q,pt 



9 2 
TyTz 



9 ? 



)fei<lz) for the transverse phonon. Here ^ = 8.5 



eV is the deformation potential; d — 5.31 g/cm^ is 
the mass density of the crystal; Vgi = 5.29 x 10'^ m/s 



J 



{vst = 2.48 X 10^ m/s) is the velocity of the longitu- 
dinal (transverse) sound wave; kq = 12.9 denotes the 
static dielectric constant and Kqo = 10.8 denotes the 
optical dielectric constant; and ei4 — 1.41 x 10^ V/m 
represents the piezoelectric constant. = 35.4 meV 

is the LO phonon frequency, and the AC phonon spec- 
tra flQ\ are given by fiq; — VgiQ for the longitudi- 
nal mode and flQt = VgtQ for the transverse modej2^ 
Nqq^^x — [exp(/3f2qg^^A) — 1]^^ rcprcscnts the Bose distri- 
bution. 



The spin-conserving electron-electron Coulomb scat- 
tering terms are given by 



dt 



dt 



= {-2n J2 Ke./(4-q - 4 + 4' - 4'-q) [(1 - /k-q,.)/k,.(l " /k' ) /k' -q..' 
q,k',(T' 

+ ^PkPk-q(/k',<T' - fk'-q.a') + ^Pk'Pk'-q(/k-q,fT - /k.tr) | - |k ^ k - q, k' ^ k' - q| , (A7) 

= {-TT Ke,/(4-q - 4 + 4' - 4'-q) (/k-q. + lPk + /k, - iPk-q) (/k' - /k'-q,a') 

q,k' ,cr' 

+ Pk[(l - ./k',ff')/k'-q.<T' - R.c(pk'Pk'_q)] - Pk-q[/k',<T'(l - /k'-q,<T') - Rc(/3k'Pk'-q)] } 

- |k ^ k - q, k' ^ k' - q| . (A8) 
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